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ABSTRACT 

In order to contribute to the general effort aiming at the improvement of our knowledge about the 
physical conditions within the Broad Line Region (BLR) of Active Galactic Nuclei (AGN), here we 
present the results achieved by our analysis of the spectral properties of a sample of 90 broad line 
emitting sources, collected at the Sloan Digital Sky Survey (SDSS) database. By focusing our attention 
mainly onto the Balmer series of hydrogen emission lines, which is the dominant feature in the optical 
wavelength range of many BLR spectra, we extracted several flux and profile measurements, which we 
related to other source properties, such as optical continuum luminosities, inferred black hole masses, 
and accretion rates. Using the Boltzmann Plot method to investigate the Balmer line flux ratios as a 
function of the line profiles, we found that broader line emitting AGN typically have larger Ha/H/3 
and smaller FLy/H/J and H<5/H/3 line ratios. With the help of some recent investigations, we model 
the structure of the BLR and we study the influence of the accretion process on the properties of the 
BLR plasma. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: emission lines — galaxies: Seyfert 
— line: profiles 



1. INTRODUCTION 

The dominant features in many Active Galactic Nu- 
clei (AGN) spectra are broad emissio n lines which origi- 
nate in the Broad Line Region (BLR) (|Osterbrock 111989k 
iKrolik 1 119991 : iPeterson I l2003h . The BLR can poten- 
tially provide a useful probe of the central source, thus 
understanding its physics and kinematics is a crucial 
step in the investigation of AGN. There are three rea- 
sons: (i) the kinematics of the BLR is probably con- 
trolled by the central source, with the competing ef- 
fects of gravity and radiation pressure; (ii) the BLR 
reprocesses the X-ray/UV energy emitted by the con- 
tinuum source, consequently the broad emission lines 
can provide indirect information about this part of the 
spectrum; (iii) there is indication that the parameters 
of the broad lines (coming from the BLR) may be re- 
lated to other fundamental properties of the source. 
Most of the recent BLR studies have been focused on 
geometries, sizes, and correlations between the BLR 
kinematical prop erties and the general characteristics of 
AGN (see e. g. ISulentic. Marziani. &: Dultzin-Havcan I 



120001 IPopovic et al. I 12004 iKaspi et al. I 120051 etc.). 
while other works have been devoted to recon- 
structing the physical conditio n s in the BLR emit- 
ting gas (Ikaspi fc Netzer 1 119991: IPopovic I [20031 120061 : 
iKorista fc Goad 112004 IVeron-Cettv et al. II2006L etc.). 

The broad line strength, width, and shape are pow- 
erful tools for gas diagnostics in th e different parts o f 
the emitting region of AGN (e. g. lOsterbrock1ll989h . 
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However, there is the problem that the broad emission 
lines are complex and that they are probably coming 
from at least two regions with different kinematical and 
physical conditions. Furthermore, the broad emission 
line profiles of some AGN m ay be explained with a two- 
component model (see e. g. IPopovic et al. [[2002, 2003, 
[2003: iTuTet al. 1120061 : iBon et al. Il200lh ~ 

Here we describe an extensive investigation of the spec- 
tral properties of the broad components in the Balmer 
emission lines, combined with the study of those other 
source physical parameters that we were able to infer 
from our sample of spectra. Exploiting some recent 
results that have been achieved by mean s of the Re- 
verberation Mapping (R M) technique (cfr. IKaspi et al. I 
120051 : iBentz et al.~ll2006lL we try to determine the central 
source masses and luminosities throughout the sample 
and to costrain their role in controlling the BLR struc- 
ture and dynamics. 

The paper is organized according to the following plan: 
in § 2 we describe the selection of our sample and we sum- 
marize the required reduction process before performing 
our measurements; in § 3 we report the extraction of our 
measurments and their related uncertainties; § 4 will list 
our results, with a discussion of our findings and their 
limits, while in § 5 our conclusions are given. 

2. SAMPLE SELECTION AND PROCESSING 

The set of spectra for our data sample (see an example 
in Fig. []}, has been collected at the spectral database 
of the 3 rd data release from the SDSS. Observations are 
performed as survey campaigns at the 2.5m //5 modi- 
fied Ritchey-Chretien altitude-azimuth telescope, located 
at the Apache Point Observatory, New Mexico (USA) 5 . 
Data are obtained with a spectrograph, whose sensor uses 
a mosaic made up by four SITe/ Tektronix 2048 x 2048 
CCDs, which covers a wavelength range running from 
3800A to 9200A. A system of 640 optical fibers, each 

5 http://www.sdss.org/dr3 
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Fig. 1. — Multiple Gaussian decomposition of the Balmer line 
profiles in the spectrum of SDSS J1203+0229. Here we denote 
with the appropriate labels all the emission lines that we took 
into account. The observed spectrum (thick continuous line) has 
been artificially shifted upwards for a better comparison with the 
narrow (short dashed), broad (long dashed), and Fe II (dotted) 
contributions. The dot-dashed line in the bottom part of each 
panel gives the fit residuals. 

having an aperture of 3", subdivides the telescope's field 
of view, so that each exposure yields 640 spectra corre- 
sponding to as many areas in the sky. The spectral res- 
olution R of the observations ranges from 1850 to 2200. 
According to the purposes of our work, we searched the 
SDSS database for sources corresponding to the following 
requirements: 

1. objects had to be located at redshift z < 0.4, in 
particular it was required that the entire profiles of 
lines belonging to the Balmer series were covered 
by the available spectral range; 

2. only spectra where the Balmer series was clearly 
recognized, at least up to the US, have been con- 
sidered; 

3. the profile of a broad component had to be de- 
tectable for each Balmer line; 

4. such profiles had not to be affected by distortions, 
due, for example, to bad pixels on the sensors, as 
well as by the presence of strong foreground or 
background sources. 

The preview spectra provided by the database retrieval 
software were manually inspected, looking for the objects 
in better agreement with our requirements, until a num- 
ber of 115 sources were chosen from approximately 600 
candidates examined in various survey areas. Subsequent 
inspection of the spectra collected within the database 
led to the rejection of 25 objects, which were affected by 
problems that could not be detected in the preview anal- 
ysis. Therefore our resulting sample includes the spectra 
of 90 variously broad line emitting AGN, corresponding 
to ~ 15% of the candidates that we examined and lo- 
cated in the range 0.024 < z < 0.368, with an average 



redshift of 0.119. A significant fraction of our sample (45 
spectra) h as been se l ected from the collection of sources 
studied bv lBorosonl (|2003l ). 

The SDSS database provides users with pre-processed 
material, therefore spectra retrieved from the survey are 
already corrected for instrumental and environmental ef- 
fects, including the sky emission subtraction and the cor- 
rection for telluric absorptions. Calibration of data in 
physical units of flux and wavelength is also performed. 
Consequently, our preliminary reduction simply needs 
to take into account a correction for Galactic Extinc- 
tion, which we estim ated using an empirical selective ex- 
tinction function (see lCardelli. Clavton. fe Mathis Ill989h 
that was computed for each spectrum on the ba- 
sis of the Galactic Extinction coe fficients given by 
ISchlegelTl^ inkbein er. fe Davis I (|1998f ) and available at 
the Nasa Extragalactic Database (NED) e , and the re- 
moval of cosmological redshift. Since our interest lies on 
the investigation of the BLR properties, we had to iden- 
tify the broad line components in spectra. In principle, 
the BLR signature can be isolated if we know which con- 
tributions are introduced by the underlying continuum 
of both the AGN and its host and by the Narrow Line 
Region (NLR). In many cases large samples are dealt 
with by means of automatic data processing techniques, 
which can be particularly useful in a statistical sense, but 
they neglect most of the peculiar properties of the single 
sources. In our case, intrinsic differences among the var- 
ious objects certainly did not simplify the task. Broad 
and narrow lines, indeed, usually have blended profiles, 
whose final shape may critically depend on many circum- 
stances, such as the orientation of our line of sight onto 
the source, the presence of absorbing material, or the 
amount of signal originating outside the source itself. For 
this reason we undertook the task of manually identify- 
ing the BLR contributions in spectra, rather than relying 
on automatically collected determinations. We used the 
IRAF software for reduction and analysis of our data. 
At first, we performed a continuum normalization, fit- 
ting the underlying continuum shape of each spectrum 
in the rest frame wavelength ranges typically running 
between 3750 - 3850A, 4200 - 4250A, 4700 - 4800A, 5075 
- 5125A, 5600 - 5700A, 6150 - 6250A, and 6800 - 7000A, 
which were not affected by significant line contamination. 
The shape of the spectral continuum was reproduced by 
means of spline functions with order ranging from 3 to 
7 in the splot task of IRAF. As shown in Fig. [IJ sub- 
traction of the narrow line components has been achieved 
by means of multiple Gaussian fits to their profiles. Us- 
ing the narrow [O III] feature at 5007A, we extracted a 
template profile for forbidden lines in each spectrum and 
the other narrow lines were required to have compatible 
widths. The line intensity ratio for the [N II] doublet at 
A A 6548, 6584 has been fixed to 1 : 3. In a certain num- 
ber of cases, before proceeding with our measurements 
of fluxes and profiles, we had to remove the contribution 
due to the Fe II multiplets. According to the su gges- 
tions given by iVeron-Cettv. Jolv. fe Veron I (|2004f ). the 
removal of spectral contributions from Fe II can be per- 
formed by scaling and smoothing a template Fe II spec- 
trum, which w as previously extrac ted from the spectrum 
of I Zwicky 1 (jBotte et al. Il2004[ ). We scaled the tem- 

6 http:/ /nedwww. ipac.caltech.edu/ 
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Fig. 2.— Spectral decomposition of SDSS J1307-0036. In this 
plot we show the spectrum of the object in the rest frame wave- 
length range used to perform the fit (thick continuous line), the 
AGN component (thin continuous line), the host galaxy contribu- 
tion (long dashed line), and the fit residuals (short dashed line). 

plate according to the Fe II features seen in spectra and 
we measured the resulting fluxes in order to estimate the 
Fe II contribution. However this step was unnecessary 
for many objects in the sample, since the Fe II emission 
lines gave negligible contributions to the fluxes or they 
were undiscernible above the average noise intensity fluc- 
tuations in these sources. 

2.1. Host galaxy correction 

Light from targets observed by the SDSS spectrographs 
is collected within sky areas of fixed aperture, which, 
in the case of AGN, will bring to various contributions 
from the host galaxy to the total flux recorded by in- 
struments, essentially depending on the object redshift 
and on the relative importance of the host with respect 
to its AGN. A possible way to account for the influence 
of host galaxies on the fo rmation of the resul t ing sp ec- 
tra is the one proposed bv lVanden Berk et al. I (|2006f ) on 
the basis of the Karhun en-Loeve Transforms described in 
IConnollv et al. I (|1995f ) . This technique assumes that the 
total spectrum of an AGN and its host galaxy may be the 
result of a sum of orthonormal components, which make 
up a set of eigenspectra, arranged in a linear combination 
such as: 

n m 

5(A) = J> • MX)] + J^lhj ■ BjM], (1) 
i=i j=i 

where we call 5(A), Ai(X), and Hj(X) the total spectrum, 
the i th AGN eigenspectrum, and the j th host galaxy 
eigenspectrum, while a% and hj are respectively the AGN 
and host galaxy eigencoefhcients for the corresponding 
components, which do not depend on wavelength. The 
possibility to estimate the host galaxy contribution in 
the resulting spectrum stems from the use of appropriate 
sets of eigenspectra, like those computed by lYip et al~l 
(|2004al lbl). both for galaxies and AGN, from large sam- 
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Fig. 3. — Relative importance of host galaxy contamination to 
the optical continuum in the spectra of our sample. The distribu- 
tion does not show any particular dependence on the object redshift 
or its luminosity. 

pies of galaxy and AGN SDSS spectra. We used an it- 
erative x 2 minimization process, which was run over the 
rest frame wavelengths covering the range from 3700A to 
7000A, to assign the appropriate eigencoefficient values 
for the decomposition of ou r spectra. In a similar analy- 
sis to the one performed bv IVanden Berk et al. I ([2006f l. 
we found that most of the variations within our sample 
could be accounted for by introducing five galaxy and six 
AGN components, because a smaller number of eigen- 
spectra left very large residuals, whilst a larger number 
took to a dangerous overfitting of noise fluctuations. The 
decomposition of our spectra usually yielded a reduced 
residual of 0.05 < % 2 < 2 in the continuum. In Fig. [2] we 
give an example of a spectral decomposition obtained by 
combining separately all the AGN and host components. 
The figure illustrates that this method is actually able to 
fit the spectral continuum in good detail, although sig- 
nificant residuals are left in the regions corresponding to 
the line cores, where the highest order variations of our 
sample are carried out. The spectral contributions iden- 
tified as being due to the host were therefore removed 
from our subsequent measurments of the AGN proper- 
ties and they have been used to estimate the relative 
importance of the AGN and its host in the emission of 
the observed luminosities. Fig. [3] illustrates the amount 
of host galaxy contaminations detected in our sample in 
the form of a histogram. When we studied the distri- 
bution of these non-AGN components, we found no clear 
evidence for a redshift or luminosity dependence. Indeed, 
taking into account objects at larger distances, there is 
a higher probability to include more powerful AGN and 
this may balance the effect of light collected from a wider 
fraction of the host. 

3. SPECTRAL ANALYSIS 

Once the broad line components had been fairly iso- 
lated in the spectra, we performed a number of measure- 
ments in order to estimate fluxes, FWHM, and FWZI. 
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The fluxes of the Balmer lines were measured several 
times and the uncertainties to be associated to the line 
flux determinations were computed as: 



I Ccont 



■ F 



F' 



(2) 



where <7 con t represents the rms due to noise fluctuations 
estimated in the continuum close to the emission lines, 
F and /50 are, respectively, the total fluxes of the lines 
and their half maximum intensities, while crp is the stan- 
dard deviation of our multiple flux determinations with 
different choices of the continuum intensity level. A list 
of flux ratios for various Balmer lines, with respect to H/3 
from different sources, is given in Table [T] Their error 
bars are obtained as: 



&Rj - Rj ■ 



. Fi 



Uh/3 



(3) 



where Rj is the ratio of j Balmer line with respect 
to H/3. Similar flux measurments have been performed 
on the [O III] emission line at 5007A and they are also 
reported in Table [TJ 

In order to estimate the FWHM and FWZI, we re- 
stricted our attention only on the strongest spectral fea- 
tures of the BLR, namely Ha and H/3, and the [O III] 
5007 line of the NLR. After choosing a zero emission in- 
tensity level, we computed the line peak intensities and 
we used the inferred values to define a half maximum 
intensity. Half widths at half the maximum and at zero 
intensity were taken both on the red and the blue line 
wings, respectively where the line profiles crossed the half 
and the zero intensity levels. Such determinations were 
repeated several times for each line, taking into account 
different guesses to the zero intensity levels, and they 
have been eventually averaged together, thus providing 
the mean values with their standard deviations. We cor- 
rected the emission line profiles for instrumental broad- 
ening, assuming that the profiles were affected according 
to: 

^obs = Ant + C- (4) 



where Wob s and W lnt are the observed and intrinsic line 
widths, while W- ms characterizes the instrumental broad- 
ening, which, in terms of velocity units, may be expressed 
as W ins = c/R w 150km s -1 . ' 

To estimate the continuum luminosities, we performed 
averaged flux measurements over the rest frame wave- 
length range running from 5075A to 5125A. After the 
correction for Galactic Extinction had been taken into 
account, the fluxes, that we derived both in the con- 
tinuum and in the selected emission lines, were con- 
verted into specific luminosities, making the assump- 
tion of isotropic emission of radiation from the sources. 
Here we used the cosmological redshift as a distance es- 
timator, in the framework of a model characterized by 
H = 75km s" 1 Mpc" 1 , n m = 0.3, and Q A = 0.7. We 
assumed the sources included in our sample to have op- 
tical luminosities given by £5100 = 5100 • L\, where L\ 
is the specific continuum luminosity measured at 5 100 A, 
and their bolometric luminosities, which we needed in 
order to guess the accretion rate onto the central source, 
to be roughly ten times as much. 



4. RESULTS 

Since we are going to consider the influence of the 
central source of an AGN onto the BLR, we must nec- 
essarily confront the problems due to its currently un- 
resolvable structure. The conversion of observational 
data into physical parameters requires some assump- 
tions to be made about the BLR structure and dy- 
namics. At present, the most trusted interpretation 
of the broad emission line profiles, which we can of- 
ten detect in AGN spectra, invokes fast orbital motions 
of a photoionized plasma in the gravitational field of 
a Super Massive Black Hole (SBH). Matter accretion 
onto the black hole provides an energetic continuum of 
radiation that interacts with the gas distribution sur- 
rounding the central source and it eventually results in 
the observed spectra. Concerning the real behaviour 
of such gas distribution, there are many aspects which 
still have to be unambiguously clarified. Recent inves- 
tigations (e. g. IVestergaard. Wilkes, fc Barthel I l2000t 
iNikolaiuk. Czernv. fc Ziolkowskv Il2005f) suggest that the 
BLR geometry should be considerably flattened, rather 
than spherical, and its motions should occur mainly in 
an orbital configuration, as opposed to radial infalls or 
outflows. 

Given that the BLR is dynamically affected by the 
gravitational field of the central source and by its 
radiation pressure, the broad line emitting medium 
must be controlled by the balance of these forces. 
As a consequence, the mass of the central engine 
and its energetic output should be tightly related 
to the size and velocity field measured in the BLR 
(IWandel. Peterson, fc Malkan I 119991: IWu et al. I 12004 



iPeterson et al. 1 120041 : iKaspi et al. 1 I2005D . In order to 
estimate the size of the BLR, a number of AGN have 
been studied with the Reverberation Mapping (RM) tech- 
nique, which exploits the time lags elapsed since the ion- 
izing continuum variations and the corresponding line 
responses to map the radial distribution of the line emit- 
ting gas. Although this method is currently one of the 
most powerful available probes of the structure of the 
BLR, it is very expensive, because it requires large mon- 
itoring campaigns for each object under study. Using 
the available data it is nonetheless possible to investi- 
gate how the BLR radii change as a function of the AGN 
luminos ities in various spe ctral ranges. In the optical 
domain IKaspi et al. I (|2005l ) argued that the size of the 
BLR scales with the source luminosity according to: 



R 



BLR 



101t-days 



= (2.23 ±0.21) 



AL A (5100A) 



10 44 erg-s" 



0.69±0.05 



(5) 



More recently (see e. g. iBentz et al. II2006D . this kind of 
relationship has been reviewed for different RM sources 
and it has been suggested that the power law index of 
Eq. (5) may be consistent with the value of a = 0.5, 
predicted on the basis of some simple photoionization 
calculations. In the following we adopt the latter inter- 
pretation to infer the source parameters from our mea- 
surments. 

Assuming the BLR to be in virial equili brium (see e. g. 
IWoo fc Urrv~ll2002t ISulentic et al. 1120061 . and references 
therein, for a discussion of implications) and to have an 
isotropic motion pattern, the latter approximation being 
forced by the lack of multi- wavelength information about 
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Fig. 4. — Line width distribution of the sample, plotted as a 
comparison of the FWHM in H/3 (top panel) and the correspond- 
ing value in [O III] A5007 (bottom panel). The thin continu- 
ous histograms show the distribution of NLSls with FWHMjr/3 < 
2000km 9 , while the thick continuous histograms represent broad 
line emitting AGN. 

the real flattening of the BLR and its incli nation with 
respect to our line of sight (see for exa mple lBianll2005l 
for a discussion of the BLR inclination), we estimate the 
black hole masses by means of: 



Mbh = ^ |fWHM h , 



(6) 



Given the optical luminosities and the inferred black hole 
masses we estimate the accretion rates onto the central 
mass with: 

Ll '" A M Edd ^, (7) 



M 



r\c" 



^Edd 



assuming that the source bolometric luminosities can be 
roughly expressed as LboI ~ lO-^sioo and that the effi- 
ciency of the accretion process is 77 = 0.1. 

In the framework of these assumptions we find that 
our sample covers a black hole mass range between 
5.59 • 10 5 M Q and 2.96 • 10 8 M Q , with the inferred ac- 
cretion rates running from 0.024 to 0.835. We give a 
summary of the results achieved by our calculations in 
Table 2. Though most of our estimates to the physical 
and structural properties of the BLR in these objects 
have been inferred by means of simplified assumptions 
and empirical relations, so that they might in principle 
be prone to large uncertainties, the sample, as a whole, 
provides a useful reference frame, where the actual prop- 
erties of realistic sources can be expected to span. Since 
the selection of our sample is limited in redshift by the 
requirement of the Balmer series falling into the avail- 
able spectral range, it lacks the most powerful broad line 
AGN. However, as we show in Fig. [H where we compare 
the FWHM of H/3 and the [O III] forbidden line, the 
sample appears to be fairly well distributed, although 
there are no objects featuring very broad lines. Among 
the 90 spectra of the sample, 20 match the classical def- 
inition of Narrow Line Seyfert 1 (NLS1) galaxies, with 



FWHMh/3 < 2000km s , while the others are spread in 
the range of 2000km s" 1 < FWHM H(3 < 7000km s -1 . 
NLSls are typically different from other Seyfert galax- 
ies under many points of view. They are commonly 
observed to have super-solar metallicities in their nu- 
clei and they often show strong Fe II multiplets, to- 
gether with a s oft slope and high variability i n the X- 
rays (see e. g. iBoller. Brandt, fe Fink I Il996f h More- 
over they are usually observed to have relatively high 
Fe / [O III] intensity ratios and strong soft X-ray ex- 
cesses. Many explanations have been proposed by var- 
ious authors to give an interpretation of the properties 
found in NLSls. Simply assuming that the line broad- 
ening depends on the gravitational potential of the SBH 
and that the BLR is not systematically larger in these 
sources, Boiler et al. (1996) suggested that NLSls are 
powered by comparatively low mass black holes, with 
respect to normal Seyfert 1 galaxies. It has also been 
propos ed that they could be young AGN in a gr owing 
phase (|Mathur, Kuraszkiewicz, and Czernv 1 120011) . that 
the narrow line profiles could be th e effect of a partially 
obscured BLR ([Smith et al~l l2002f) or that they could 
arise fro m an extremely flat geomet ry seen at low incli- 
nations (|Osterbrock fe Pogge I [r985), but the polariza- 
tion signatures which we should expect in the last two 
cases have not been reported so far. Many of the char- 
acteristic properties of NLSls could not be further in- 
vestigated in our sample, because of its limited spectral 
coverage and of the usually not ideal s/n ratios, which, 
in many cases, did not allow us to actually measure the 
strength of Fe II emission. In conclusion, we found that 
NLSls in our sample, located at an average redshift of 
z = 0.064, have on average the Lpeii/isioo = 0.029 ± 
0.013, L Fo ii/i [0 iii] = 4.30 ± 3.08, and L FeU /L H p = 
16.82 ±7.81 luminosity ratios, while the other AGN have 
ipeii/isioo = 0.024 ± 0.020, L FeI i/i [0 ni] = 4.36 ± 5.31, 
and LFeii/^H/3 = 13.09 ± 11.66. In the following we kept 
the distinction of NLSls from the other AGN. 

4.1. The continuum source and Balmer Line properties 

At present it is strongly believed that the main inter- 
actions between the central power source of an AGN and 
the surrounding emission line regions occurs via radiative 
processes. A direct measurement of the AGN continuum 
luminosity on the spectra, however, is not straightfor- 
ward because of the several uncertainties about the ac- 
tual AGN Spectral Energy Distributions (SED), which 
can even be affected in various manners b y stellar con- 
tribu tions from the AGN host galaxy (see iBentz et al. 1 
l2006l for instance) . On the other hand, some clues about 
the reliability of the optical continuum intensity as an es- 
timator of the ionizing radiation field strength come from 
the observation of how the emission line intensities, which 
should be controlled by the amount of ionizing radiation, 
correlate with the estimated continuum luminosities. In 
Fig.O we show the distribution of our sample in the case 
of [O III] A5007 as a function of the optical source lumi- 
nosities. The behaviour that we find is described by the 
relation: 

log(i[o in]) - (0.820±0.052).log(L 5100 ) + (5.710±2.270), 

(8) 

with a correlation coefficient R = 0.712 and a probabil- 
ity to occur by chance P < 10~ 6 . Most of the scatter 
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Fig. 5. — Correlation between the [O III] A5007 line luminos- 
ity and the optical continuum luminosity measured in the spectra. 
NLSls have been identified with filled circles, while the other broad 
line emitting sources are represented by open circles. The cross lo- 
cated in the upper left corner of the diagram shows the median 
uncertainty of our measurements. The straight line fit is achieved 
by means of Eq. (8). Intrinsic absorption and starlight contribu- 
tions might be responsible for the scatter. 

observed in this relation can be accounted for in terms of 
intrinsic source absorptions or stellar light contributions 
which could have been not perfectly removed. 

Looking at the Balmer line intensity ratios over the 
sample, where we found the median values of Ha/H/3 = 
3.45 ± 0.65, H 7 /H/3 = 0.45 ± 0.08, and fW/H/3 = 
0.25 ± 0.08, we compared the properties of the sources 
with their associated line profiles. The plots shown in 
Fig. [6] illustrate the Balmer line flux ratios as a function 
of the FWHM measured in the spectra. The Balmer 
line flux ratio measurements are prone to quite large un- 
certainties, presumably reflected by the resulting scat- 
ter, but the plots suggest the existence of a weak rela- 
tion between the Balmer decrement and the width of line 
profiles. Although the correlations are not particularly 
strong (having R = 0.226, P = 3.34 • 10~ 2 in the case of 
Ha/H/3, while R = -0.482, P = 1.71 • 10~ 6 for H 7 /H/3, 
and R = -0.510, P < 10~ 6 for H5/H/3), they apparently 
strengthen with increasing line widths, but this may be 
simply due to the small number of such objects. Almost 
no particular correlations are, instead, observed among 
the Balmer line flux ratios themselves, with the only ex- 
ception of the ratios of H5/H/3 and ITy/H/3, whose values 
are related by: 

log(iWF H/3 ) = (1.253 ± 0.133) • log{F u ^/F u0 )- 

-(0.171 ± 0.047), (9) 

with a correlation strength of R — 0.709 and P < 10~ 6 . 
These fin dings agree very well with earlier results ob- 
tained by Rafanclli | (|1985f) . who, while investigating the 
relationships among the Balmer line intensities and pro- 
files in a sample of 12 Seyfert 1 galaxies, recorded a cor- 
relation like the one given in Eq. (9) (with a slope of 
1.24) and the absence of similar relations involving the 
other emission lines. 



Fig. 6. — Panels from the upper left to the lower left show 
the Balmer line flux ratios with respect to H/3 as a function of 
the FWHM measured in spectra. We observe an averagely larger 
Balmer decrement in objects with broader line profiles. The bot- 
tom right panel illustrates the only correlation observed between 
H5/H/3 and H7/H/3. We plot with open circles the broad line 
emitting sources and with filled ones NLS1 galaxies. The median 
uncertainty estimates are given by the crosses shown in each plot. 

4.2. Global Baldwin Effect 
The occurrence of a power law slope like L[ Q mi oc 
£5100° with a = 0.820 ± 0.52 in Eq. (8) ind icates the 
presence of a global Bald win Effect in this line (|Baldwin I 
ll977l:IKong etaDl2006h . As we look for the situation of 
the broad Balmer line components, however, we find lit- 
tle or no evidence for global Baldwin Effect, since the 
observed emission line luminosities are related to the op- 
tical continuum through: 

log(-LHa) = (0.989 ± 0.031) • log(L 5 i 00 ) - (0.784 ± 1.346) 

(10a) 

log(L H/3 ) = (1.029 ± 0.030) • log(L 5 i 00 ) - (3.025 ± 1.289) 

(106) 

log(i H7 ) = (0.976 ± 0.029) • log(L 5 ioo) - (1-063 ± 1.226) 

(10c) 

log(L H i) = (0.931 ± 0.030) • log(L 5 ioo) + (0.631 ± 1.299), 

(lOd) 

with correlation coefficients of R = 0.937, 0.967, 
0.974, and 0.972, spanning from Eq. (10a) to Eq. (lOd), 
and P < 10~ 6 in all cases. The situation of the 
Balmer lines is summarized in Fig. [7] Here we note 
that, although there are works which actually detected 
a significant intrinsic Baldwin Effect, particularly in 
the case of H/3, in long tim e covering observations of 
some well studied objects ([Gilbert fc Peterson I 120031 : 
iGoad. Korista. &; Knigge 1 12004( 1 , no strong indications 
have been found in literature concerning the existence 
of a glob al Baldwin Effect in the broad Balmer lines 
Csee e. g. lYeel 119801: iBinette. Fosburv. fc Parker~lll993l : 
lOsmer fc Shieldsl fl99flh . 

4.3. Balmer Decrements and Boltzmann Plots 

The broad line flux ratios estimated in our sample, 
which we presented in Fig. [6l are a potential probe to 
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Fig. 7. — Correlations between the optical continuum luminosity 
and the H<5 (upper left panel), H7 (upper right panel), H/3 (lower 
left panel), and Ha emission line luminosities (lower right panel). 
Open circles represent broad line emitting objects, while filled cir- 
cles are NLSls. The straight lines are the best fits of Eq. (10) and 
the error bars show the median uncertainties associated with our 
measurements. 

study how the physical properties of the BLR influence 
the Balmer decrement. To investigate this effect we need 
to choose a specific parameter that should be used as an 
estimator of the Balmer decrement and of its shape in 
each object. A suitable possibil ity is to apply the Boltz- 
mann Plot (BP) method (see iPopovic I l2003t 120061 for 
description) to the broad line components of the Balmer 
series. Introducing a normalized line intensity with re- 
spect to the atomic constants involved in the transition: 

-'n = a ' V-L-U 

g u ■ A u i 

where F u i is the measured flux, while X u i, g u , and A u i 
are the line wavelength, the upper level statistical weight, 
and the spontaneous radiative decay coefficient respec- 
tively, it can be shown that the line intensities of a spe- 
cific series of transitions in an optically thin plasma de- 
pend on the excited level's energy according to: 

4 (12 > 

in which we denoted with Nq the number density of 
the radiating species, Z the partition function, E u the 
excited level energy, and T e the excitation temperature, 
while £ is the radial extension of the emitting region, and 
h, fcs, and c are the fundamental constants of Planck, 
Boltzmann, and the speed of light. If the high excitation 
stages (with n > 2) in hydrogen are well described by the 
Saha-Boltzmann distribution, Eq. (12) takes the form of: 

log/ n =B- A-E u , (13) 

where A is called the Temperature Parameter because, 
in the situation depicted above, it would be A = 
loge/(fc B T c ) « 5060/T e [K -1 ]. 

When we apply the BP to the Balmer series of a BLR 
spectrum, the resulting slope A clearly depends on the 



Fig. 8. — Four examples of Boltzmann Plots applied to the 
spectra of objects belonging to the sample: in the upper left panel 
no straight line fit is achieved (class iv); in the upper right panel 
only a poor fit is performed (class iii); the bottom left panel gives 
a good fit, but He could not be detected in the spectrum (class 
ii); finally the bottom right panel shows a straight line fit to the 
normalized intensities of the Balmer series up to He (class i). 

Balmer decrement, though there is a substantial differ- 
ence between these two properties, since, for instance, 
we find objects with a rising BP slope (A < 0), although 
this does not imply Ha to be really fainter than the other 
Balmer lines. The shape of the Balmer decrement, too, 
influences the BP, which can be dramatically limited by 
the restrictions involved in the assumptions about the 
physical conditions within the plasma. Intrinsic redden- 
ing effects of the AGN environment, too, may influence 
the method, enhancing low order lines with respect to 
the high order ones, thus taking to a steeper straight 
line fit. Therefore, in presence of intrinsic reddening, the 
observed BP temperatures would be probably lower than 
the actual values, although there are indications that the 
straight line fit to the normalized intensities is not qual- 
itatively affected (|Popovic | [2003). According to whether 
the observed spectra could be described by a relation like 
Eq. (13), we studied the results of BP and we found the 
occurence of four most common situations, correspond- 
ing to (i) a good fit of the Balmer series up to the He emis- 
sion line, (ii) a good fit where, however, the He line could 
not be detected, (iii) a poor fit to the observed lines, and 
(iv) no possible straight line fit of the line normalized 
intensities. Examples for each one of these classes can 
be found in Fig. [8] We report the BP slopes and an indi- 
cation of the proper class in Table 1, while, in Fig. [9j we 
show the distribution of the inferred BP slopes as a func- 
tion of the BLR velocity fields. Here we see that the BP 
slope increases on average w ith the broad line wid ths, as 
it was previously reported (jPopovic 1 120031 . [2006) . The 
total degree of correlation is quite weak (R = 0.381, 
P = 2.27 • 10~ 4 , in the case of the relation between A 
and FWHM, while it is R = 0.470, P = 3.28 • 10~ 6 for 
that between A and FWZI), but it grows significantly if 
only the cases where a reasonable BP fit is achieved are 
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Fig. 9. — BP slope A as a function of the BLR velocity fields 
inferred by both the FWHM and FWZI of the Balmer lines. Here 
we plot with open circles the sources where the straight line fit 
covered the series up to He (i) and with open squares those objects 
having a fairly good fit up to H<5 (ii). Filled triangles are spectra 
with a poor fit (iii) and finally asterisks are spectra where the 
Balmer line normalized intensities could not be suitably described 
by a straight line (iv). The crosses in the lower right corners of 
both panels are the median uncertainties of measurements. 

taken into account (R = 0.453, P = 1.08 • 10" 3 for A and 
FWHM, R = 0.579, P = 1.30 • 10~ 5 for A and FWZI). 

4.4. Balmer decrements and Eddington ratios 

The rate of the accretion process that produces the 
ionizing continuum of radiation in AGN is a crucial pa- 
rameter affecting the physical conditions of the emis- 
sion line regions. Many theoretical works suggest that 
the SED of an ionizing radiation field depends on the 
rate and the efficiency of the accr e tion process (e. 
iNetzer. Laor. fc Gondhalekar |[l99l iKong et"ai~ll200C 
which, in turn, affects the ioniza tion status of the BLR 
and its stability (|Nicastro | |2000). Unfortunately the de- 
terminations of accretion rate, estimated from black hole 
mass and luminosity, may be considerably affected by 
model dependent assumptions. In our case, the situation 
is shown in Fig. 1 101 where we compare with our sam- 
ple the expected distribution of AGN powered by black 
holes with masses in the range 10 4 M Q < Mbh < 10 9 M Q 
and accreting at the labelled values of their Eddington 
ratios. We note that this model does not predict any 
particular difference in the accretion rate of NLSls and 
of the other sources. The accretion rates that we infer 
by assuming a stru ctural model computed according to 
iBentz et al. I (|2006[ ) show a rather small degree of corre- 
lation (R = -0.314, P = 2.17- 10~ 3 ) with the strength 
of the Balmer decrement, as summarized in the BP slope 
A. As we illustrate in Fig. [Til this means that the high 
order lines appear to be stronger when the accretion rate 
is large, with, perhaps, a flattening where the line ratios 
approach the expected atomic values. In order to under- 
stand the relationship between the accretion rate and A, 
we should clarify the role of this parameter. As men- 
tioned above, in the appropriate circumstances the tem- 




km s 



Fig. 10. — Bolometric luminosities vs. FWHMjj/3 in the sam- 
ple. Filled circles represent NLSls, while open circles are broad 
line AGN, and the cross in the upper left corner gives the median 
uncertainty of the measurements. The curves track the expected 
locations of black holes with masses between 1O 4 M0 (bottom left 
side of the diagram) and 10 9 Mq (upper right region) accreting at 
the labelled values of M/M^ddi if their emission were isotropic, the 
ac cretion efficiencies w ere r) = 0.1, and the BLR scaled as described 
in IBentz et~al~l pOOlf) . 

perature parameter A is related to the thermodynam- 
ics of the emission line region. Although the physical 
conditions of the BLR are such that a thermodynamic 
interpretation could not be generally straightforward, it 
might be possible to explore some fraction of the BLR 
with this technique. Especially in the objects where we 
found fairly good fits to the BP, the dependence of A 
on the line profiles and Eddington ratios could be a con- 
sequence of the influence of ionizing radiation onto the 
emission line plasma. In this case the occurrence of large 
temperature parameters in broad line emitting sources 
is an indication of low ionization, maybe because of a 
shielding of the broad Balmer line emitting gas, while the 
anti-correlation between the Balmer decrement and the 
Eddington ratio might be a consequence of high dust col- 
umn densities in our sightline towards the objects with 
low accretion rate. It is possible that the temperature 
parameter A tracks the influence of the ionizing radia- 
tion field onto the broad line emitting plasma, therefore 
its relations with the line profiles and Eddington ratios 
would suggest a stronger ionization in narrow line and in 
high accretion rate sources, but a more detailed expla- 
nation of this result would only be possible by a deeper 
understanding of its role as a thermodynamic diagnostic 
tool. 

5. DISCUSSION AND CONCLUSIONS 

The main purpose of this work is to collect a large sam- 
ple of flux and profile measurements of the broad emis- 
sion line components in the Balmer series and to compare 
the results with those AGN properties that we were able 
to infer from the available data. The paper is intended 
as a starting point in a larger investigation aiming at the 
identifiaction of possible hints about the physical condi- 
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FlG. 11. — Balmer line flux ratios of H8 (upper left panel), H*y 
(upper right), and Ha (lower left) together with the BP slope pa- 
rameter (lower right panel) as a function of the Eddington ratio, 
with the same symbols as in Fig. [9] There appears to be a weak 
anti-correlation between the accretion rate and the Balmer decre- 
ment, with the high order lines being stronger in objects working 
at high Eddington ratios, resulting in modest BP slopes. 

tions within the BLR of many AGN, through the analysis 
of various broad lines. 

While we spent a great effort in performing as care- 
fully as possible our direct measurements on the broad 
emission lines, the size of our sample and the limited 
spectral coverage forced us to calculate the other source 
parameters by means of some simplified assumptions and 
empirical relations, which, in some cases, are still matter 
of debate. Using the BP method, we explored the rela- 
tion of the temperature parameter A and the broad line 
profiles. The values inferred in those objects where the 
technique achieved the best results would correspond to 
a temperature range running from 8 • 10 3 K to 21 • 10 3 K, 
with cold gas being usually associated with broad line 
emitters. A better explanation of this result, which could 
be of fundamental importance in our understanding of 
the BLR structure, requires a deep investigation of the 
properties of the temperature parameter A as a thermo- 
dynamic tool (?). The possibility of a relation with the 
inferred accretion rates, too, deserves further investiga- 
tion, since it may inform us about the different condi- 
tions of gas ionization or dust processing of the radiation 
within our point of view towards the BLR. 

Taking into account the results of our measurements 
and the limits of the related assumptions, here we come 
to the following conclusions: 

• the Balmer line flux ratios show a weak degree of 
correlation with the line profile width, in the sense 
that more pronounced Balmer decrements are ob- 
served in broader line emitting objects; 



• the optical continuum luminosity of AGN is related 
to the emission line intensity of [O III], where we 
confirm the presence of a global Baldwin Effect, 
and of the broad components of the Balmer series, 
which, instead, do not show any strong evidence of 
global Baldwin Effect; 



the shape of the Balmer decrement, studied with 
the Boltzmann Plot technique, is related to the line 
profile widths and we observe that larger slopes and 
better straight line fits can be achieved in objects 
with broader lines, where a BLR component being 
described by the assumptions of this method may 
be present; 



• we find a weak degree of anti-correlation between 
the Balmer decrement and the Eddington ratios 
measured in our sample, that could be the result 
of high dust column densities along our line of sight 
to the objects characterized by low accretion rates. 
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TABLE 1 

Line Flux and Boltzmann Plot analysis 
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± 





07 





.47 


± 





.07 


3 


,88 


± 


0. 


32 


3564 ± 131 


0. 


334 


± 





,037 


ii 


SDSSJ1659- 


-6202 


1440 ± 101 


ND 





.25 


± 


0. 


04 





.41 


± 





.05 


3 


10 


± 


0. 


24 


149 ± 27 





251 


± 





,019 


ii 


SDSSJ1717+5815 


1906 ± 144 


ND 





20 


± 


0. 


.04 





.40 


± 





.07 


2. 


,27 


± 


0, 


22 


521 ± 49 





096 


± 





,065 


iv 


SDSSJ1719+5937 


5646 ± 297 


ND 





.13 


± 





02 





.37 


± 





.04 


4 


,36 


± 


0. 


34 


1955 ± 144 


0. 


525 


± 





072 


ii 


SDSSJ1720+5540 


3722 ± 253 


ND 





,24 


± 


0. 


04 





40 


± 





.04 


3 


,54 


± 


0. 


26 


5491 ± 213 


0. 


323 


± 





,016 


ii 


SDSSJ2058-0650 


1116 ± 191 


ND 


0. 


.20 


± 





.06 





39 


± 





.10 


2 


,57 


± 


0. 


48 


296 ± 30 


0. 


204 


± 





036 


iii 


SDSSJ2349-0036 


1837 ± 176 


ND 





.34 


± 





.06 





,54 


± 





.08 


3 


.05 


± 


0. 


33 


519 ± 45 


0. 


146 


± 





,043 


iv 


SDSSJ2351-0109 


522 ± 52 


ND 


0. 


.26 


± 





.04 





.47 


± 





.09 


3 


,06 


± 


0. 


35 


186 ± 17 


0. 


201 


± 





,013 


iii 



a Fluxes arc given in units of 10 17 crg cm 2 s' 
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TABLE 2 

Line profile and source property estimates 



Object Name 


FWHM a 


FWZI a 


-^5100 b 




M BH 


d 




M/M 


Edd 


srtss mm 3-nQ5i 


1 1 99 


4- 


1UU 


c:ne:i _i_ 1 coo 
uyuj. zlz J-Uoz 


q Qfi 4_ 1 7 

O . OU ZLZ U.Li 


11.35 


± 
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8.08 ± 0.93 
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4_ 
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± 


62.92 
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± 
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± 
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Q QQ 4- 08 

o.oy zlz u.uo 


11.40 


± 
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± 
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U 
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± 
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± 
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940 
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4- 
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qdqq toi 1 94-0003 
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4_ 


1 7Q 
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± 


1.19 
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± 
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U 
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ZL 
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± 
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± 
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TABLE 2 

Line profile and source property estimates 
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a FWHM and FWZI arc measured in km s 1 £5100 arc gi vcn in units of 10 42 crg s 1 C figL^ arc in units of lO^cm^ M BH arc expressed in units of IO^Mq 



